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Quantum effects at motion of ultracold neutrons in a magnetic 
trap are actively studied now. In particular, "magnetic tube" 
traps [I] are used for investigation of such effects. According to 
PP, the periodic set of magnetic traps can be developed. Traps 
in the form of "magnetic tube" are used for focusing of atomic 
beams and described in details in [2j on the basis of the classic 
motion equations for a particle, possessing the magnetic moment, 
in a magnetic field. But quantum effects, which appear at neu- 
tron motion in the magnetic trap, can not be described by classic 
motion equations. The typical radius of the magnetic trap, which 
is being designed, is p ~ 10~ 4 cm pQ|. Thus, due to large wave- 
length (A ~ 1CT 5 -7- 10~ 6 cm), motion of the ultracold neutron in 
the magnetic trap in the direction, transversal to the axis of the 
magnetic tube, appears quantized. 
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Quantum effects were studied for neutron motion in a gravita- 
tional trap [3]. Neutron motion in such a trap is one-dimensional. 

But in contrast to [3] , transversal motion of neutrons in the " mag- 
netic tube" trap with respect to it axes is two-dimensional. As 
noted below, neutron motion in such a trap is determined by spin- 
independent effective potential energy, which is proportional to 
the magnetic field strength square. 

Thus, let us consider motion of the ultracold neutrons in the tube- 
shaped magnetic trap. The Schrodinger equation, which describes 
the ultracold neutron in such a trap, is as follows: 

l~^r ~ ^B(r}\ tt(f) = E*(r) (1) 
where M is the neutron mass, p is the neutron magnetic moment, 
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r is the coordinate of the neutron, B is the magnetic induction 
in the point r , a = (a x , a y , a z ) are the Pauli matrices. 

Suppose the axes z to be directed along the tube axes. Neutron 
motion along the tube axes is quasi-classic. Tube length (dimen- 
sion in direction z) significantly exceeds its radius. Therefore, 
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B{f) inside the tube can be considered as independent on z, i.e 
B(f) = B(x,y) = B(p,(p), where p and <p are the radius and 
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azimuth angle in the cylindrical coordinate frame. 

The equation (CO) in the cylindrical coordinate frame looks like as 
follows: 
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Particle motion along z is free and is described by the quantum 
number k z . Therefore, the wave function of a neutron in the tube 
can be expressed as: 

oo 

\&(?) = <f>(p)e ik * z = £ <£> m (p)e imip e ik * z (3) 

m=— oo 

Therefore, to find the spectrum of neutron states in the tube the 
following equation should be solved: 
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$m + ^E^„ m '(p)$m'(p) = 0(4) 
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where x 2 = k 2 -k 2 , k 2 = B mml {p) = ± J e i( - m '- m ^B(p, <p)d<p 

Let us find the spectrum of neutron states in the trap in the first 
order of perturbation theory. In this case the diagonal matrix 
element of the interaction energy, i.e. the average value of the 
magnetic field strength B in the trap, determines the correction, 
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which defines spectrum. But the averaged value of B in such a 
trap is zero. Therefore, the nonzero contribution to the spectrum 



of the neutron energy levels in the tube trap appears in the second 
order of the perturbation theory. 

For further analysis it is convenient to transform differential equa- 
tion (J2J) to the homogeneous integral equation: 

i£(r) = J G Q {r^)V{^^dy (5) 

where Go(r, f*) = ~2^2 V_p| ' ^ s the Green function of the equa- 
tion (ED and the potential energy V = — /io\B(r). Substitution of 
the expression for ^(r) (0) to the right part of equation (0) gives: 

m{f) = J G (r, r , )V{f , )G {r' } r")V(f")y(f")d 3 r'd 3 r" (6) 
This equation is equivalent to the following expression 
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A r - E\ #(f) + J V{r)G {r, 7)V{7)^{7)d 3 r' = 0(7) 



Motion of the neutron along z is quasiclassic, therefore the qua- 
siclassic expression for the Green function can be used: 

G (f, = -i-j^W ~ P)°(z ~ z)e ik{z '~ z) (8) 



where 6(z' — z) is the Heaviside unit function (6(z' — z) = 1, when 
z' > z and 9(z' — z) = 0, when z' < z). 
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With the help of (ED the expression (J7D converts to 

^A r - e\ *(r) - ^V(P)V(P) Je- ikz e ik *'V(p, z')dz' 
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It should be mentioned that V(p)V(p) = /i 2 (aB)(aB) = fi 2 B 2 (p). 



With the help (ED the equation (J2j) can be expressed as 



h " A p - \ *W + V eff (p)$(P) = (10) 
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where 



v Mfi 2 B 2 (p) 



is the effective potential energy. 



According to ( 1111 ) neutron motion in the magnetic tube trap, 



which is formed by the alternating magnetic field, is determined 
by the squared magnetic field strength and does not depend on 
the neutron spin direction. 

The equation ( TTUft allows to find the eigenfunctions and the spec- 
trum of eigenvalues of neutron energy in the trap. 

In the cylindrical coordinate frame the equation ( TTUT) reads as 



follows: 

+ -lT- ^r + ^) <Mp) - -Zr Urn\V eff (p)\m')<S> m ,(p) = 0(12) 
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dp 2 pdp p 2 J h 2 
heie<S>(p) = E m $m(p)e^, (m\V eff (p)\m') = ± f 2 * V eff (p, ip)e^ m - m >dtp. 

Let us keep in the sum in ( fT2l ) only the term with m! = m, which 
describes the effective potential energy V e ff(p,(p) averaged over 
the azimuth angle (p. (More detailed analysis of neutrons in the 
trap could require consideration of additional terms in the sum 
with m! 7^ m.) Then ( fT2l ) looks like as follows: 
d 2 Id M 2 \ 2M 

^2 + -Q- p - ~j + ^ 2 J $m(p) ~ ^{m\V eff {p)\m)<$> m {p) = 0(13) 

According to [?] magnetic field in the tube trap linearly decreases 
to the trap center, i.e. B(p) ~ p. Therefore, B 2 (p) ~ p 2 and 
(m\V e ff(p)\m) ~ p 2 . Therefore, equation ( TT31 ) can be rewritten 
as: 

rl 2 \ r) M 2 \ 

jjf + - -r + " V$m(p) = o (14) 

where ^4 is determined by 

, 2 2M. , . 2Mp 2 (m\B 2 (p)\m) . , 

The eigenvalues and eigenfunctions of the obtained equation ( IT41 ) 



are the those of the cylindrically symmetric oscillator [3j. Intro- 
ducing new variable £ = \fAp 2 one can get from ( HM the following 
equation: 

^ + £ +( -! + *-?>— (16) 
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where = -^7=7. 

Solution of this equation is as follows [?] : 

MO = (17) 

where W(£) is the confluent hypergeometric function 

W(() = F L(J3 - ^Ltl), |m| + uj . (18) 

To make the wavefunction $ TO finite everywhere, the parameter 
/? — l m ^ +1 should be integer nonnegative number, i.e. 

772 H - 1 

/? = rip, where n p > 0. (19) 

This condition ( 1191 ) determined spectrum of the transversal mo- 
tion for the neutron in the magnetic trap. 

To avoid misunderstanding, I'd like to emphasize that I consider 
spectrum of the energy levels, which are lower then the height of 



the effective potential barrier. 
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Substitution j3 = in the explicit form to ( TT9T ) provides to get: 
x 2 I ml + 1 



4v^4 
therefore 



= rip, where n p > 0, (20) 



^ / \ h 2 K 2 (n ,m) h 2 , 1—, |m| + l N , n , 
Etrans(n p , m) = v p ' 7 = 4Vi n„ + — 21 



Let us consider the particular example. The averaged field can 
be expressed as 

(m\B 2 (p)\m) = ^A p 2 (22) 
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where p max is the radius of the magnetic tube, B^ ax is the square 
averaged magnetic field at p = p max . In this case the coefficient 
A is 

A = 2M 2 p 2 (B 2 max ) 
h 2 k(k + k z )p 2 nax 

The frequency of the transverse oscillations of the neutron in the 
ground state in the magnetic tube 



ntmns = ^ = ^2 pSSS. (24) 
M HJk(k + k z )p max 



For the ultracold neutrons with A ~ 10 6 -j- 10 5 cm in the tube 
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with p max ~ 10 4 cm the product ^Jk(k + k z )p max « fc/? moa! > 
1. Therefore, for the magnetic field \f{B^J) = 10 4 gauss the 
frequency O trans ~ 10 4 10 5 s . Suppose that during the time 
r the neutron moves with the speed 10 2 cm/s in the tube of 1 
cm length (time of neutron flight in the trap is about r ~ 10~ 2 
s), then the parameter fl trans T ^> 1, i.e. neutron executes many 
oscillations in the trap. 

Thus, neutron motion in the magnetic tube is described by the 
energy of the transversal motion E tra ns an d the energy of longi- 
tudinal motion Ei ong = -^f. The sum of these energies is equal 
to the neutron energy before it flies into the trap E = 

E = Ef rans + E{ OTl g. (25) 



The equation ( |25| ) can be used to find k\\ for the neutron in the 



certain state of transversal motion. 

When the energy of transversal motion is comparable or exceeds 
the height of the effective potential barrier, the above expressions 
for the wavefunctions and E trans cannot be applied. 

The typical angle of neutron incidence on the magnetic tube trap 
providing the energy of transversal motion to be less than the 



height of the effective potential barrier can be expressed as: 



where is the frequency of Larmor precession of neutron spin in 



the field v(B^ aa .). The expression (J26J) was obtained considering 
k m k z (because k z converges to k with energy growth). 

According to ( 125T) the neutron moving with the speed 1 m/s in 
the trap with the maximal field B « 10 4 gauss can be captured 
to the trap with the potential V e ff at the angle d < d CT « 1 rad. 
When neutron speed is w = 10 m/s, the angle $ cr « 10 -2 rad. 
Energy growth results in reducing of # cr ~ ~ ^ 

Neutron passing through the magnetic tube trap is similar to the 
process of axial channelling of a charged particle in a crystal. 
Study of the angular distribution for neutrons passing through 
the trap provides to find spectrum of the neutron levels in the 
trap. 

Actually, experiments for neutrons(atoms) passing through the 
magnetic tube trap can be considered as experiments studying 
neutron (atom) channelling in such a trap. When the traps are 
periodically distributed in the plane (x, y) the analogy becomes 
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ever more evident. 

It should be also mentioned that bent magnetic tube traps can 
be used for focusing wide beams of ultra cold neutrons (similar 
focusing of particles channelled in a crystal). 
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